The Active Traveling Wave in the Cochlea 
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A sound stimulus entering the inner ear excites a deformation of the basilar membrane which 
travels along the cochlea towards the apex. It is well established that this wave-like disturbance is 
amplified by an active system. Recently, it has been proposed that the active system consists of a 
set of self-tuned critical oscillators which automatically operate at an oscillatory instability. Here, 
we show how the concepts of a traveling wave and of self-tuned critical oscillators can be combined 
to describe the nonlinear wave in the cochlea. 

PACS numbers: 87.10.+e, 47.35.-K, 05.45.-a 



The sounds that we hear are detected within the cochlea 
of the inner ear, a fluid-filled duct which is coiled like the 
chamber of a snail's shell. This compartment is parti- 
tioned along almost its entire length by the basilar mem- 
brane (BM) . Incoming sound waves set the BM into mo- 
tion and its minute vibrations are monitored by special- 
ized sensory hair cells 0. The pioneering experiments 
of von Bekesy Q, which were conducted on cadavers, 
demonstrated that sound excites a traveling wave on the 
BM, whose amplitude reaches a peak at a place which de- 
pends on the frequency. This suggested that the cochlea 
acts as a spatial frequency analyzer. When it became fea- 
sible to measure the BM response of living specimens, a 
marked difference from von Bekesy's results was revealed. 
The sharpness of filtering was greatly enhanced and the 
response displayed pronounced nonlinear behavior close 
to resonance 0, E| • These observations, together 
with the discovery that ears spontaneously emit sounds 
at specific frequencies Q , provided direct evidence for an 
active nonlinear amplifier in the cochlea Q, as had been 
foreseen by Gold |£g. But just how the ear could reliably 
employ an active process without suffering from unstable 
behavior has long been a matter of concern. 

An active amplifier which overcomes this difficulty has 
recently been outlined 0, 0, Q| . It has been proposed 
that the cochlea contains a set of dynamical systems, 
each of which is maintained at the threshold of an os- 
cillatory instability by a self-tuning mechanism. Poised 
at this critical point, on the verge of vibrating, each sys- 
tem is especially responsive to periodic stimuli at its own 
characteristic frequency. The concept of self-tuned crit- 
ical oscillators 10] can account for the main features of 
hearing: sharp frequency selectivity, extreme sensitivity 
and wide dynamic range; and also for interference effects 
such as two-tone suppression and the generation of com- 
bination tones 12]. In this letter, we marry the concept 
of critical oscillators with the physics of the traveling 
wave to provide a unifying description of active cochlear 
mechanics. 

Cochlear waves. — The basic physics of cochlear waves 



may be described most succinctly by a one-dimensional 
model [HQ EE HE!- The BM separates the cochlear 
duct into two channels which are connected at the apex 
by a small aperture, the helicotrema. A sound stimulus 
impinging on the oval window, at the base of the cochlea, 
causes changes in the pressures P\(x,i) and P2(x,t) in 
both channels. Here t is the time and x is the position 
along the cochlea, with the oval window at x — and 
the helicotrema at x = L. The pressure gradients induce 
longitudinal currents J\{x,t) and J2{x,i), which flow in 
opposite directions in the two channels. We define the 
relative current j = J\ — Ji and the pressure difference 
p = Pi — Pi- Then the balance of pressure gradients and 
inertial forces in the fluid may be written 

pdtj = -bld x p , (1) 

where p is the fluid mass density, I is the height of each 
channel, and b is the breadth of the BM. The conservation 
of fluid volume implies that a variation in the current 
along the cochlea must be accommodated by a movement 
of the cochlear partition. We describe such deformations 
of the BM by its height h(x, t) as a function of time and 
position. Then the conservation law is 

2bd t h + d x j = . (2) 

Combining this with Eq. we obtain an equation for 
the BM acceleration 

2 P bd 2 t h = d x [bld x p] . (3) 

The pressure difference p acts to deform the BM. If the 
response is passive (eg. in the dead cochlea), the response 
relation close to the basal end, where the stiffness K(x) 
of the BM is high, takes the simple form 

p(x,t) = K(x)h(x,t) , (4) 

for small disturbances. Eqs. © & 0) together yield 
a linear wave equation for the pressure, with local wave 
propagation velocity 
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Critical oscillators. — In the active cochlea, the passive 
response is amplified by a force-generating system. This 
system comprises a set of mechanical oscillators which 
are supported on the BM, and which are positioned in 
such a way that they can drive its motion. The char- 
acteristic frequency cu r (x) of the oscillators is a function 
of position along the membrane. In general, such oscil- 
lators could either vibrate spontaneously and thus gen- 
erate motion in the absence of a stimulus, or they could 
be quiescent and behave like a passive system. A par- 
ticularly interesting case arises at the boundary of these 
two regimes, when every oscillator operates exactly at 
the critical point where it undergoes an oscillatory in- 
stability. Automatic regulation to this critical point - 
or Hopf bifurcation — can in general be achieved by us- 
ing a robust self-tuning mechanism based on local feed- 
back control [10(. If the BM contains such self-tuned 
critical oscillators, its deformation h in response to pres- 
sure differences across the membrane p has characteristic 
properties as a function of frequency and amplitude, and 
nonlinear amplification occurs. 

In order to describe this system, we first consider an 
individual oscillator. Its characteristic response to peri- 
odic forcing at frequency u> can be written in a general 
form as [li| 

p = A(uj)h + B\h\ 2 h . (6) 

Here, h and p are the Fourier amplitudes at the forcing 
frequency and A and B are complex coefficients. This 
expression follows from a systematic expansion in the os- 
cillation amplitude h which is valid close to the Hopf bi- 
furcation (comparable to a Landau expansion of the free 
energy of thermodynamic systems near a critical point). 
Proximity to an oscillatory instability thus automatically 
provides for generic nonlinearities. The dominant nonlin- 
earity is cubic, a result that follows from time-translation 
invariance. The linear response coefficient A vanishes at 
the characteristic frequency u> r of the oscillator so that, 
at this particular frequency, the response becomes purely 
nonlinear for small amplitudes. 

Thus if we focus on a particular location x of the BM, 
its response displays a nonlinear resonance when the fre- 
quency of the stimulus is equal to the local characteris- 
tic frequency u> r (x) of the oscillators. The shape of the 
resonance, for nearby frequencies, can be described by 
expanding the function A(u>) in powers of u — u> r (x). For 
frequencies that differ substantially from the local char- 
acteristic frequency, on the other hand, we expect the 
active system to contribute little to the BM response. In 
particular, when u> = 0, the BM deflection is governed 
only by its passive stiffness, according to Eq. J2J. We 
now assert that the simple functional form 



A(x,w) — a(u> r (x) — uj) 



(7) 



where a is a real constant, captures the essential features 
of this BM response. Clearly it satisfies the requirement 



that the linear response coefficient at location x can be 
expanded about uj r (x). Secondly, it indicates that the 
passive stiffness is proportional to the characteristic fre- 
quency: K(x) = A(x,0) = aoj r {x). This corresponds 
well with experimental data. The frequency-place map 
and the elasticity of the BM have been carefully mea- 
sured. Characteristic frequency and stiffness both de- 
crease approximately exponentially with distance along 
the cochlea, falling by about two orders of magnitude 
from base to apex We therefore supplement Eqs. 
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ui r (x) = loqc 
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to obtain the full position-dependent response of the BM. 
We take the coefficient B, describing the nonlinearity 
close to resonance, to be a purely imaginary constant, 
B = i[3. This simple choice ensures that Eq. © has no 
spontaneously oscillating solution for p = 0, as required 
at the critical point. 

Active traveling waves. — Combining Eq. J3J) for the BM 
acceleration with the response of an active membrane, 
described by Eq. 10, we obtain a nonlinear wave equation 
for the BM deformation. In frequency representation, it 
reads 



-2pbu> 2 h = d x 



bid. 



(A(x,L 



)h + B\h\ 2 h 



(9) 



The complex solutions of this equation h(x) — H(x)e' l ^ x ^ 
describe the amplitude H and the phase cj> of the BM dis- 
placement elicited by a periodic stimulus with incoming 
sound pressure p(x = 0, t) = p(0)e luJt . 

For small pressures, the nonlinearity described by the 
cubic term in Eq. (0 is significant only close to the res- 
onant place which, inverting Eq. is x r = dln(u>o/Lo). 
Far from this characteristic place, we obtain a linear wave 
equation which can be solved in the WKB approximation 
|14l ll(il| . The approximate solution can be expressed as 



h(x) 



(Cj r (x) -W) 3 / 4 
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i I dx' q(x') 
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with local wave vector 



>\ 1/2 



laj (ui r (x)-uj) 



1/2 



(10) 



(11) 



At the basal end of the cochlea, x < x r , q is real and the 
solution is a traveling wave with a position-dependent 
wave vector. As the wave propagates, its wavelength di- 
minishes and its amplitude builds up, until it approaches 
the place of resonance. In the immediate vicinity of the 
characteristic place, A decreases according to Eq. (J7J) 
while h increases. Thus the cubic term in Eq. rapidly 
becomes more important than the linear term. This cuts 
off the divergence in Eq. (|10fl and leads to a strongly 
nonlinear BM response. The wave peaks at x = x p < x r , 
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where the response displays the characteristic nonlinear- 
ity of critical oscillators, h(x p ) ~ p(x p ) 1 ^ 3 0,0. From 
Eq. (JTDJl we find that p(x p ) ~ A{x p ) 1 / A p(Qi), while the 
crossover from linear to nonlinear response implies that 
A{x p ) ~ |S||ft-(a;p)| 2 . We thus find that the peak ampli- 
tude has a power law response 

h(x p )~ P (oy (12) 

as a function of the stimulus pressure at the base, with an 
exponent v = 0.4. At positions beyond the characteristic 
place, x > x r , the wave vector q becomes imaginary, 
indicating the breakdown of wave propagation. The BM 
displacement decays very sharply in this regime. 

Numerical solutions. — Full solutions to the nonlinear 
wave equation, Eq. Q, can be obtained numerically. It 
is most convenient to solve the equation for the pressure 
p(x) which satisfies 

- t ip=(A + Bu(p))d 2 J , (13) 

with /i = 2pto 2 jl and where we have assumed for sim- 
plicity that b and / are constant along the cochlea. The 
variable u — H 2 is the squared deformation amplitude 
and is a nonlinear function of p. Indeed, it follows from 
Eq. JHJ that u(p) is the unique real and positive root of 
the cubic equation 

\p\ 2 = \A\ 2 u+ (A*B + AB*)u 2 + |B|V . (14) 

Eq. (|13|) for p can be solved, starting from x = L and 
integrating towards x = 0. As a boundary condition, 
we impose zero pressure difference at the helicotrema, 
p(L) = 0, because the two cochlear channels are con- 
nected there. A second boundary condition specifies the 
value of d x p at x = L. By varying this pressure gradient 
at the helicotrema. we find solutions that correspond to 
waves entering the cochlea at x = with different pres- 
sure amplitudes p(0). The profile of BM displacements 
can then be obtained from the solution p{x) via 



Basilar membrane response. — Examples of traveling 
waves are displayed in Fig. 1 for two different sound lev- 
els and varying stimulus frequencies. Waves initiated at 
x = propagate with growing amplitude and decreasing 
wavelength until they reach a point of resonance, beyond 
which they decay rapidly. At 40 dB SPL, the resonance 
is sharp and the peak response occurs at a location very 
close to the characteristic place x = x r , where the fre- 
quency of the active oscillators is equal to the stimulus 
frequency. At 80 dB SPL, the resonance is much broader 
and the location x ~ x p of maximal response shifts to- 
wards the base, in agreement with experimental observa- 
tions @|. 

The response at a particular location on the BM ex- 
hibits the qualitative properties of cochlear tuning that 
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FIG. 1: The active traveling wave on the BM for three differ- 
ent frequencies (/ = 370 Hz, 1.3 kHz & 4.6 kHz) whose corre- 
sponding characteristic places are x T /L — 0.25 (top), 0.5 (cen- 
ter) & 0.75 (bottom). Instantaneous BM displacement h(x,t) 
(black lines) and wave amplitude H (x) (gray lines) are shown 
for two stimulus amplitudes, characterized by sound pres- 
sure level (SPL). At 40 dB SPL, a sharp resonance occurs at 
the characteristic place, where the critical oscillators actively 
drive the response. The resonance is broader at 80 dB SPL, 
and the peak shifts towards the base of the cochlea. Param- 
eter values: Cochlear dimensions L — 35 mm and I = 1 mm; 
fluid density p — 10 3 kg/m 3 ; parameters governing the fre- 
quencies of the active oscillators u>o = 10 s s _1 and d = 7mm, 
providing a frequency range of 100Hz-16kHz. We choose 
a — 10 4 Pa/ms, which implies a volumetric stiffness of the 
BM varying in the range 6 x 10 6 - 10 9 Pa/m. There is only 
one free parameter in our calculations, f3 = 10 23 Pa/m 3 , which 
we choose to fit the nonlinearity of the response according to 
sound pressure level (SPL). For simplicity, it is assumed that 
the middle ear raises sound pressures by 20 dB, independent 
of frequency. 



have been observed experimentally |3J, U, |£j . Fig. 2a dis- 
plays the gain w|ft,|/|p(0)| of BM velocity, obtained from 
our numerical solutions, as a function of stimulus fre- 
quency for different sound levels. At low frequencies the 
response is linear and the gain is independent of the stim- 
ulus amplitude. As the stimulus frequency approaches 
the resonant frequency, the response becomes nonlinear 
and the gain diverges as the SPL declines. At higher 
frequencies, the response drops precipitously. The mag- 
nitude of the BM displacement, shown in Fig. 2b, is typ- 
ically several nanometers at resonance, in quantitative 
agreement with experimental data Q. The phase <j> of 
the traveling wave at a particular location on the BM is 
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lem with this interpretation is that, in order to obtain the 
observed range of characteristic frequencies, the stiffness 
of the BM would have to vary by more than four orders of 
magnitude from base to apex. The measured variation is 
only a factor of one hundred 0, [2^ • This difficulty is cir- 
cumvented by our theory, where the range of frequencies 
at which the BM resonates is determined by the frequen- 
cies of the oscillators that are ranged along it, and is not 
governed by the stiffness or the inertia. Some models of 
the active cochlea are very specific and rely on additional 
degrees of freedom, secondary resonances or time-delayed 
feedbacks |2l|. Such descriptions lack the simplicity and 
generality of our approach and miss t he g eneric nature 
of the power-law nonlinearities [lfl fill \2?fy conferred by 
the Hopf bifurcation. 

In this letter it has been our aim to provide a concise, 
coherent interpretation of a wide variety of observations, 
rather than a detailed fit of individual data. Consider- 
ing that our model incorporates only one free parame- 
ter whose value is not determined by independent mea- 
surement, the qualitative agreement with a diverse set 
of experimental data is striking. We have not sought to 
specify the physical nature of the active oscillators. The 
electromotility of outer hair cells has been implicated in 
active amplification in the mammalian cochlea [H Eil . 
but the motile response of hair bundles may also play a 
role. Indeed, the hair bundles of frog hair cells have re- 
cently been demonstrated to behave as Hopf oscillators 
[261 127| . Because the response of self-tuned critical oscil- 
lators is generic, our analysis remains valid whatever the 
physical basis of force generation. 

We thank M. Magnasco, P. Martin, E. Olson, J. Prost 
and G. Zweig for stimulating discussions. T.D. is a Royal 
Society University Research Fellow. 



FIG. 2: Response of the BM at a fixed position as a func- 
tion of stimulus frequency, for different sound pressures. The 
characteristic frequency of the active oscillators at that lo- 
cation is 5 kHz. (a) Velocity gain, i.e. BM velocity divided 
by sound pressure of the stimulus, (b) Corresponding BM 
displacement, (c) Phase difference between stimulus and BM 
oscillation. 



displayed in Fig. 2c. As observed experimentally, it de- 
creases with increasing frequency — gradually at first, 
but more abruptly as resonance is approached — and 
then varies only little at frequencies higher than the char- 
acteristic frequency. 

Discussion. — In many recent cochlear models, sections 
of the BM are considered to behave as inertial oscilla- 
tors which are either lightly damped (in passive models) 
14Lll5Lll6| or driven by internal forces (in active models) 
lit l2(i |21| . The characteristic frequency at a particular 
location is then the local resonant frequency of the BM, 
which varies as the square root of the stiffness. A prob- 
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